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Orbital    Transfer  with  Minimum  Fuel 

•    School 

V.  E.  BLEICK  and  F.  D.  FAULKNER* 
U.  S.  Naval  Postgraduate  School,  Monterey,  Calif. 

A  note  in  this  Journal,  Ref.  1,  discussed  the  problem  of  scheduling 
the  direction  p  of  constant  momentum  thrust  of  a  rocket,  which  loses 
mass  at  a  constant  rate,  so  that  it  transfers  to  a  known  earth  satel-  • 
lite  orbit  in  minimum  time  T  after  launching.  A  numerical  solution 
was  obtained,  using  rectangular  coordinates,  for  the  case  of  fixed 
launching  conditions.  The  method  of  Ref.  1  is  extended  here  to  solve 
the  problem  of  orbital  transfer  of  such  a  rocket  with  minimum  fuel 
consumption.  All  of  the  symbols,  units,  and  end  conditions  of  Ref.  1 
are  used  here  without  redefinition. 
Statement  of  the  Problem 

The  time  of  flight  T  in  minimum  fuel  transfer  must  be  longer  than 
in  the  minimum  time  transfer  of  Ref.  1,  unless  these  two  trajectories 
turn  out  to  be  identical.  This  implies  at  least  one  interruption  in 
rocket  thrust  during  minimum  fuel  transfer.  The  problem  solved  here 
assumes  exactly  one  such  interruption,  i.e.  launch  at  t=0,  thrust 
interruption  at  t=tq  ,  thrust  resumption  at  t=t~  ,  and  final  thrust 
termination  at  transfer  t=T.  The  problem  of  minimum  fuel  transfer  is 
equivalent  to  the  Lagrange  calculus  of  variations  problem  of  requir- 
ing the  integral 


J  =  /n  (f+Xcp-j+nq^+TT^+pcp^dt  (1) 


-X 

to  be  stationary,  where  f  is  the  fuel  consumption  rate,  X,,n,rr,p  are 
continua  of  Lagrangian  multipliers,  and  cp^  =cp_=cp^=cp  =0  are  the  first 
order  equations  of  rocket  motion  of  Ref.  1 .  The  f  function  and  the 


Received  September    ,  1963.  This  work  was  supported  by  the  Office 
of  Naval  Research.  Aid  of  the  Computer  Facility,  U.  S.  Naval  Post- 
graduate School,  is  acknowledged. 

*  Professors,  Department  of  Mathematics  and  Mechanics. 


rocket  thrust  per  unit  remaining  mass  function  a  are  defined  as  fol- 
lows: For  0<t<t1  ,  f=1  and  a=cM/(l -Mt) g.  For  t1<t<t2  ,  f=0  and  a=0. 
For  t2<t<T,  f=1  and  a=cfi/[l  -M(  t+t.,  -t2)  ]g.  Note  that  for  t2<t<T 
3a/dt-  =  -da/dt2  =  ga  /c.  The  varied  time  subinterval  end  points  in 
Eq.O)  are  taken  as  t- +At-  ,  t2+At2  and  T+AT.  The  vanishing  first 
variation  6J  and  its  partial  integration  are  computed  as  in  Ref .  1 . 
The  coefficients  of  6u,6v,6x,6y ,6p  in  6J=0  give  the  Euler  Eqs.(2) 
and  (3) ,  consisting  of  the  adjoint  equations 

X+TT=0  ,  £  +  p=0  , 

(2) 
*+g!  x^g2-xH=0    ,  p+gly\+g2yM=0    , 

and  the  control  equation 

tan  p  =  |i/X   .  (3) 

The  coefficient  of  AT  in  6j=0  gives,  with  the  aid  of  Eq.(3),  the 

transversality  condition 

(a.A)T  =  (aA)T  =  1  (4) 

2   2  1  /2 
where  the  adjoint  vector  A  =  iX+jn,  A=|A|=(X+ii)/  f  and 

a  =  a(i  cosp  +  j  sinp) .  The  coefficients  of  At-  and  At~  in  6j=0  give, 

with  the  aid  of  Eq.(4),  the  corner  conditions 

H(t.,)  =  [aA]J  -  JA   a2Adt  =  0,    H(t2)  =  0.  (5) 

Eqs.(5)  are  equivalent,  by  the  definition  of  a,  to 

Mt.,)  -  A(t2)  (6) 

and,  by  partial  integration,  to 

H 

L      aAdt  =  0  .  (7) 

Numerical  Solution 

Let  X  .  ,u .  ,rr .  ,p  .  ,  i=1,2,3,4,  be  four  linearly  independent  solutions 
of  the  adjoint  Eqs.(2)  corresponding  to  the  columns  of  the  matrix  E(t) 
of  Ref.  1.  The  control  angle  p  of  Eq.(3)  is  defined  by 

tanp  =  (n1+l^2+mH3+^4)/(k1+lk2+mX3+nX)-  (8) 

-  2  - 


and  its  variation  6p  is  obtained  in  terms  of  6l,6m,6n  by  total  dif- 
ferentiation as  in  Ref .  1 .  The  Bliss  fundamental  formulas  are  obtained 
by  assuming  that  a  solution  of  the  rocket  motion  equations  cp^=cp  =cp_=cp  =0 
has  been  found,  corresponding  to  Eq.(8),  which  does  not  necessarily 
satisfy  the  terminal  conditions  at  t=T  or  the  corner  condition  Eqs.(5). 
Using  this  solution  and  holding  T  fixed,  but  allowing  t.,  and  t2  to 

vary,  find  the  variation  of  the  vanishing  matrix  integral 

"T 


J0    [cp1  ,cp2,cp3,cp4]E(t)dt  =  0  (9) 


with  the  terminal  constraints  at  t=T  removed.  Since  the  columns  of 

E(t)  satisfy  the  adjoint  Eqs.(2),  one  obtains  the  system  of  Bliss 

formulas  in  the  1x4  matrix  equation 

[6u,6v,6x,6y]TE(T)  +  [G(t1 )-(apP)T]At1 

-  [G(t2)-(apF)T]At2  =  [0,6l,6m,6n]A   (10) 

where  the  matrix  A  has  been  defined  in  Ref.  1 ,  and  where  the  matrix 

G(t)  =  (apP)^  -  */   a2pFcLt  (11) 

where  the  2x4  matrix  F(t)  is  the  first  two  rows  of  E(t) ,  and  where 

the  matrix  p=[cosp, sinp] .  Substitution  of 

[6u,6v,6x,6y]T  =  [U-u,V-v,X-x,Y-y]T  +  [U-u,V-v,X-x,Y-y]TAT     (12) 

into  Eq.(lO)  gives  four  of  the  required  six  Newton-Raphson  equations 

for  the  determination  of  AT , At- , At2 ,6l ,6m,6n  on  the  varied  trajectory. 

The  remaining  two  equations  attempt  to  satisfy  the  corner  condition 

Eqs.(5)  on  the  varied  trajectory.  Involved  here  are  the  differentials 

da  =  6a   +  adt 

9  (13) 

=    (ea/dt1)At1    +   (5a/dt2)At2   +   (gaVc)dt 

and  dA   =   6A    +  Adt 

(14) 
as   {.0, 61,6111,611  JF^'    +    (Xcosp+|Jsinp)dt 

where  the  primes  on  F  and  p  indicate  matrix  transposition.  Use  of 

Eqs.(6)  and  (14)  yields  the  Newton-Raphson  equation 

A(t^Atq  -,A(t2)At2  -  [0,6l,6m,6n][F'p']t2  =  A(t2)  -  a(^)  .    (15) 

-  3  - 


Use  of  Eqs.(l3)  and  (14),  and  the  first  of  Eqs.(5),  yields  the 
Newton-Raphson  equation 

(aA)TAT  +  (K-aA)t  At1  -  K(t2)At2  +  [0,6l  ^m^njG1  (t1 )  =  -HU.,)   (16) 
where 

K(t)  =  f[a2/N]^  -  2(f)2^Ta3Adt  .  (17) 

The  iteration  to  successive  varied  trajectories,  using  Eqs.(10), 
(12),  (15)  and  (16),  may  be  carried  out  as  in  Ref.  1.  Two  devices 
were  used  to  stabilize  the  course  of  the  iteration.  The  first  was  to 
adjust  the  m  and  n  values  of  the  new  T,t.,  ,t2  »i  ,m,n  sextuple,  found  by 
solving  the  Newton-Raphson  equations,  to  satisfy  the  corner  condition 
Eqs.(5)  before  proceeding  with  the  next  iteration.  The  second  device 
was  to  modify  the  |_U, V,U, V9X,Y _u  terms  in  the  Newton-Raphson  equations, 
before  solving  these  equations,  so  as  to  minimize  the  sum  of  the 
squares  of  the  elements  of  [U-u, V-v^X-XyY-y]™. 

The  numerical  example  of  minimum  fuel  transfer  given  here  involves 
the  same  launching  conditions,  mass  loss  parameters  and  circular 
orbit  used  in  the  minimum  time  transfer  of  Ref.  1 .  The  results  for 
minimum  fuel  transfer  are  T=0. 353977,  t.,  =0.21 0293,  t~=0. 275349 , 
l=-0. 820196,  m=-0. 708727,  n=-1. 181390,  and  the  transfer  sector  angle 
B=0.1 89345  rad.  Since  the  minimum  time  trajectory  of  Ref.  1  gave 
T=0. 289869,  the  net  fuel  saving  in  minimum  fuel  transfer  over  minimum 
time  transfer  is  measured  by  0. 289869-0. 353977+t2~t1  =  0.000948,  or 
an  unspectacular  one-third  per  cent.  Figure  1  shows  the  trajectories 
and  thrust  directions  for  minimum  time  and  minimum  fuel  transfer. 

The  semilogarithmic  plots  of  Fig.  2  show  the  different  behavior  of 
^'versus  time  in   tb»  two  problems.  For  some  reason  thera  is  much  more 
difference  than  we  expect.  The  curve  increases  monotonically  for  min- 
imum time  transfer.  The  curve  for  minimum  fuel  shows  a  rather  char- 

U;      U       «&    _4_     .     i     : 


acteristic  shape.  It  is  large  initially  and  decreasing;  then  it 
increases,  and  then  decreases.  If  the  final  decreasing  interval 
does  not  occur,  larger  values  of  T  lead  to  lower  values  of  fuel 
consumption,  as  may  be  partially  inferred  from  Eq.  (7). 

Reference 

1 

Bleick,  V.  E. ,  "Orbital  transfer  in  minimum  time,"  AIAA 

Journal  1,  1229-1231  (1963). 

Figures 

Fig.  1   Trajectories  and  thrust  directions 

Fig.  2   A  versus  time 
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-YVAR 

=    -YV 

N(K>     =    C 

(K)     =    YV 

LAM)        15 

MU)  15 

=    90.0 


*AAA 
AA(3 
AA(3 
AA(4 
A*CL 
Y(30 

M+1) 
(M)» 
5 
P*DY 

ARS( 
STE 
K-l) 
S(  1) 
(2) 
(5) 
(4) 

5)  + 

6)  + 
XLAM 
LAM 
LAM 
AM 


(2)* 

)*AA 

)  *AA 

)«AA 

AM 

) 

*cos 

XA«X 

(J) 

J)    + 
P 
+    S 


2)    -    SINP*YC(4*M+1 ) 

AAA(M+1 J/CLAM 

AAA(M+1  )/CLAM 
A(3  J/CLAM 
A(4)/CLAM 
A (4) /CLAM 


P  +  YC(4*M+2)*SINP 
A 


(  AK( 1 , J)+2.*AK(2, J)  +  2.*AM3 ,J)+AK(4,J) 
TEP 


)/6 


EL*YVARS(9)  +  EM«YVARS(13)  +  EN*YVARS(17) 
EL*YVARS(10)  +  EM*YVARS(14)  +  EN*YVARS(18) 
**2  +  XMU«*2) 


8 

1) 


(4*M+1)*C0SP  +  YVARS(4*M+2)*SINP 

7,4 

S  (  M+  2  ) 

ARS( M+3) 

OVERG*XA*CLAM 

ARS( 30) 

9,152,159 

7,  155,153 
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79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

1  10 

1  1  1 

1  12 

1  13 

1  14 

1  15 

1  16 

117 

1  18 

1  19 

120 

121 

122 

123 

124 

125 

126 

127 

123 

129 

130 

131 

132 

133 

13u 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

157 

158 

159 

160 

161 

162 

163 

164 


1  5U 

155 

156 

157 

158 

159 

160 

161 

162 

165 

16fc 

165 

166 

167 

168 

169 

170 

171 

172 

173 

17U 

175 

176 

177 

178 

179 

ISC 

181 

182 

182 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

19U 

195 

196 

197 

198 

199 

20C 

201 

202 

203 

20U 

205 

206 

207 

208 

209 

210 

211 

212 

215 

hu 

215 
216 
217 
218 
219 
220 
221 
222 
223 
22-U 
225 
22  6 
227 
228 
229 
230 
231 
252 
235 
234 


GO  TO  165 
P(K)  =  0.0 
GO  TO  165 
P(K)  =  -90.0 


GO  TO  165 
P(K)  =  57 
IF  (XLAM) 
IF  (Xf-'U) 


2957 
161 
164, 


P(K)  =  P(K) 
GO  TO  165 
P(K)  =  P(K)  - 
XLAMDCT  =  0Y(5 
XMUOOT   =  0Y(6 
IF  (Nl+l-K)   1 
00  169 
C  I  T  1  (  M  ) 
A2LAMT1 


*  ATANF(XMU/XL 
165, 165 
162,  162 


AM) 


+  130.0 


CLAMDTl  = 
CGAL0T1  = 
CLAMT1 
C2ACGT1  = 
C3ACGT1  = 
C4ACGT1  = 
ATI  =  XA 
CGALMT1  = 
0A2LMT2  = 
CA3LMT2  = 

R 


180.0 

)    *     EL*DY{9)    + 

)    +     EL*OY( 10)     + 

85,  168,  185 
M=l,4 
=    CC(M) 
=    XA*XA*GLAM 

=    COSP*XLAMCCT    +     SINP*XMUDOT 
-    C0VERG*XA*CLAMDT1 
=    CLAM 

=  C0VERG*XA*CC(2) 
=  C0VERG*XA*CC(3) 
=  CCVLRG*XA*CC(4) 


EM*DY(13)  +  EN*DY{ 17) 
EM*DY(14)  +  EN*DY(18) 


COVERG*XA*CLAM 

YVARS(3C) 

2.0*YVARS(31 J/COVER 
DO     182       M  =  l,4 
0CIA2T2(M)     =    YVA^S( 31+M) 
IF     (M2+1-K)        190,  184,  190 
DO    185       M=1,4 
CIT2(M)     =    CC(M) 
A2LAMT2    =    AT1*AT1*CLAM 
CGALAM(K)     =    CCVERG*AT1*CLAM 
CLAMDT2     =    CCSP*XLAMDOT    +    SINP 
CLAMT2    =    CLAM 
IF     (N4-K)       206,191,206 
CLAMDT    =    CCSP*XL4MD0T    +    SINP* 
CGALDT    =    COVERG*X A*CLAMDT 
C2ACGT    =    C0VERG*XA*CC(2) 
C3ACGT    =    CCVERG*XA*CC(3) 
CUACGT    =    COVERG*XA*CC(U) 
CGALMT    =    COVERC*XA*CLAM 
QA2LMT    =    YVARS( 50 ) 
CA3LMT    =    2.0*YVARS(51 )/COVERG 
DO    20C    M=l  ,4 
QCIA2KM)     =    YVARS  (51+M) 
XR    =     SQRTF( YVARS( 5)**2    +    YVAR 
=    -YVARS ( 3) /XR**3    +    XA* 
=    -YVARS(U) /XR**3    +    XA* 
=       YVARS(I) 
YVARS(2) 


*XMUDOT 
XMUDOT 


S(4)**2) 

COSP 
SINP 


N15X2HBC15X2HT1 1 3X2HT21 3X1 HT  ) 

T 


CY(  1  ) 
DY(2) 
DY03) 
DY(4)  = 
CONTINUE 
PRINT  203 

FORMAT ( 1H06X2HEL15X2HEM15X2HE 
PRINT  21C,  EL, EM, EN,BB, Tl ,T2, 
FORMAT  (7F15.9) 
PRINT  212 

FORMAT  (  1HC6X2HN1  1 3X2 HN2 1 3X2HN5 1 3X 1 HU 1 4 X 1 H V 1 4X 1 HXl  4X1HY) 
PRINT  214,  Nl  ,N2, N5,XU(N4 ) ,XV(N4)  ,XX(N4) ,XY(N4) 
FORMAT  ( 2115, 4F15 .7) 
PRINT  216 

FORMAT  (  1H05X4HCAPU1 1 X4HCAPV1 1 X4HC APX 11 X4HCAP Y } 
PRINT  218,  (CAPV(M),  M=l,4) 
FORMAT  (4F15.7) 
1=1  ,4 


DO  22^ 
A(  I  ,1  ) 
B(  I  )  = 
DO  22U 
A( I  ,  1  ) 
B(  I  ) 
DO 
A(  I 


A(  I 

DO 

AA( 

AA( 

AA( 

AA( 

AA( 

DO 


227 
,2) 


=  0. 
CO 

J=l 

=  A< 
=  B( 
1  =  1 


,4 
I, 
I  ) 

,4 


1)  +  YVARS(4*I+J 
+  YVARS(4»I+J 


,3) 
23C 

U  J 

2,  J 
3,2)  = 
3,4)  = 
4,4)  = 
23c 


=    AT1*CIT1 ( I) 
=-ATl*CIT2 ( I) 
J  =  2,4 
)=    YVARS( 19+J) 
YVARS(22+J) 
YVARS(27 ) 
YVARSJ28 ) 
YVARS{ 29 ) 
1=1,4 


(0CIA2 
(QCIA2 


)*(DY( J)-CAPVD( J) ) 
)*(CAPV( J)-YVARS( J) ) 

T( I )-GCIA2T2(I ) )/COVERG 
T{  I)-QCIA2T2(I) )/COVERG 


)  = 


165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
173 
179 
180 
181 
182 
133 
184 
135 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 

20  8 
209 
210 

21  1 
21 
21 
21 
21 
21 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
235 
254 
255 
226 
257 
258 
239 
240 
241 
242 
243 
244 
245 
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'235  DO  236   J  =  l  ,  I  2U6 

236  AA( I , J)  =  AA( J, I)  2U7 

237  CO  239   1=1, 4  2U8 

238  DO  239   J  =  l,3  249 

239  A(  I, J  +  3)  =  AA(  I,J  +  1  )  250 

240  A(5, 1  )  =  -CGALDT  251 

241  A(5,2)  =  A2LAMT1  +  CGAL0T1  +  CA3LMT  -  CA3LMT2  -  XA*XA*CLAM  252 

242  A(5,3)  =-A2LAMT2  -(CA3LMT  -  CA3LMT2)  +  XA*XA*CLAK                   253 

243  A(5,4)  =  0CIA2K2)  -  CCIA2T2(2)  +  C2ACGT1  -  C2ACGT  254 
2UU  A(5,5)  =  GCIA2TI3)  -  0CIA2T2(3)  +  C3ACGT1  -  C3ACGT  255 
21+5  A(5,6)  =  CCIA2TU)  -  QCIA2T2(4)  +  C4ACCT1  -  C4ACGT  256 
2^6    B(5)  =  CGALMT  -  CGALMT1  -  GA2LMT  +  QA2LMT2                        257 

247  A(6, 1 )  =  0.0  258 

248  A(6,2)  =  CLAKDTl  259 

249  A{6,3)  =-CLAMDT2  260 

250  00  251   J=2,n  261 

251  A(6,J+2)  =  CITl(J)  -  CIT2(J)  262 

252  R<6)  =  CLAMT2  -  CLAMT1  263 

272  DO  27U   1=1 ,N4  283 

273  CGALAM(I)  =  CGALAM(N4)  -  CGALAM(I)  284 

274  A2LAT-HI)  =  A2LAM(N4)  -  A2LAM(I)  285 
IF  (KK)   320,320,300  263.1 

300  KK  =  KK-1  263.2 
GET  =  A{ 5,5)*A(6, 6 )-A { 5, 6  )  *A ( 6  ,5 )  263.3 
EM  =  EM  +  ( B( 5)*A (6,6)-B(6J»A(5,6) )/OET  263.4 
EN  =  EN  +  (B(6)*A(5,5)-B(5)*A{6,5) )/DET  263.5 
PRINT  301  263.6 

301  FORMAT  ( 1 H04X 3HTA U8X6HCAP LAM7X6HCGALAM8X5HA2L AM/ )  263.7 
PRINT  302,  (TAU(I ) ,CAPLAM( I) ,CGALAM(I ) ,A2LAM( I ),   1=1, N4)  263.90 

302  FORMAT  (4F13.9)  263.92 
GO  TO  37  263.95 

320  Dbl  =  ( XU(N4)-GAP V( 1 ) )/CAPV(2) 
DB2  =-( XV(N4)-GAP V(2) )/CAPV( 1 ) 
DB3  =  (XX(N4)-CAPV( 3) )/CAPV(4) 
DB4  =-(XY(N4)-CAPV(4) )/CAPV( 3) 
BB  =  BB  +  (DBl+DB2+DB3+DB4)/4. 
SINB  =  SINF(BB) 

GCSB  =  GCSF(RB) 

CAPV( 1 )  =  V*C0S8 

CAPV(2)  =-V*SINB 

GAPV( 1)    =  R*SIN8 

CAPV( 4)  =  R*CCSB 

CAPVD( 1 )=-VSGDR*S INB 

CAPVD(2)=-VSGDR*C0SB 

CAPVD(3)=  CAPV(  1) 

GAPVD{4)=  CAPV(2) 

PRINT  208  263.99 

PRINT  21C,   EL,EM,EN,BB,T1 ,T2,T  263991 

321  PRINT  21c 

322  PRINT  213,  (CAPV(M),  M=l,4) 

323  DO  ^328   1  =  1  ,4 

324  A(  I  ,  1  )  =  CO 

325  B ( I )  =  O.C 

326  DO  328   J=l,4 

327  A(I,1)  =  A(I,1)  +  YVARS(4*I  +  J)*(DY( J)-CAPVD(  J)  ) 

328  EMI)  =    B{I)  +  YVARS(4*I+J ) * ( CAPV ( J ) -YVARS ( J) ) 

253  GALL  GAUSS3  (6,  0.1E-09,  A,  AI ,  KER)  264 

254  PRINT  255,  KER  265 

255  FORMAT  (5H0KER=I1)  266 

256  IF  (KER-2)   258,257,258  267 

257  STOP  257  268 

258  DO  261   1=1,6  269 

259  DEL ( I )  =  0.0  270 
26C  DO  261   J=l,6  271 

261  DEL(I)  =  DEL(I)  +  AI(I,J)*B(J)  272 

262  DO  263  1=1 ,6  273 
265  TVAR(I)  =  TVAR(I)  +  DEL(I)  .                         274 

264  BB  =  BB  +  V*DEL(1)/R  •   275 

265  PRINT  208  276 

266  PRINT  210,  EL T EM, EN , BB , Tl , T2 , T  277 

267  IF  (T2-T)   303,303,307  278 

303  IF  (Tl-T)   304,304,307  278.1 

304  IF  (T2)   307,305,305  279 

305  IF  (Tl)   307,306,306  279.1 

306  IF  (T1-T2)   271,271,307  280 

307  GO  TO  275  280.  1* 
271  CONTINUE  282 

275  PRINT  276  286 

276  FORMAT ( 1H04X3HTAU 1 0X2 HXU 1 1X2HXV1 1X2HXX1 1 X2HXY9X6HG A PL AM7X6HC G ALAM8287 
1X5HA2LAM9X1HP/)         ■  288 

278  PRINT  279,  ( TALK  I  ) , XU { I ) , XV(  I )  ,  XX  {  I )  , XY ( I ) ,CAPLAM(  I ) , CGALAMt I ) ,    289 
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1  A2LAN (  I  ) ,P( I  )  , 
9    FORMAT  (3F13.9, 
■  STOP  280 
il  END 

SURROITINE  GAUSS3 
DIMENSICN  A(6,6) 
CO  1  1  =  1  ,N 

co  l  j=i,n 

1  X(  I, J  )=C.C 
CC  2  K=1,N 

2  X(K,K)=1  .0 

C  DO  34  L=1,N 
KP  =  C 
Z  =  0.0 

DO  12  K  =  L,N 
IF(Z-ARSF(A(K,L)) 
Z=APSF(A(K,L)  ) 
K  P  =  K 

CONTINUE 

IF(L-KP) 13,20,20 
DC  14  J=L,N 
Z=A(L, J ) 

)=A(KP 

J)  =  Z 
J=l  ,N 

,J) 

)=X(KP 


1=1 ,N4) 
13.2) 


(N, EP,A,X,KER) 
X(6,6) 


290 
291 
292 
293 
294 
295 


)  11,  12, 12 


J) 


J) 


)  )-EP)50,'50,30 

34 


36,32 


A(L,  J 

A(KP, 

CC  15 

Z  =  X(L 

X(L,  J 

X(KP,  J)=Z 

IF(ABSF( A(L,L 

IF(  L-M31  ,34, 

IP1=L+1 

CC  36  K=LP1,N 

IF( A( K,L ) )32, 

RATIO  =  A(K,L)/A(L,  L) 

DO  33  J=LP1,N 

A(K,J)=A(K,J  J -RAT  IO*A(L, 

DO  3  5  J=l ,N 

X(K,J)=X(K,J)-RATIC*X(L, 

CONTINUE 

CONTINUE 

DO  4  3  I=1,N 

II=N+1-I 

CC  43  J=1,N 

S  =  0.0 

IF(  II-N)41 

IIP1=II+1 

DO  42  K=  I  I 

S=S  +  A{  I  I  ,K 

X(  II,  J)  =  (X 

KER=1 

RETURN 

KER  =  2 

END 

END 


J) 
J) 


,43,43 

P1,N 

)»X (K,J  ) 

{  II, J)-S)/A(  II  , 


II) 


296 


0004C 
00050 
00C6C 
0C070 
0C08C 
0009C 
0C10C 
00110 
00120 
001  3C 
0C1U0 
00150 
00160 
00170 

coiao 

00190 
0020C 
00210 
00220 
0023C 
002^0 
00250 
00260 
C0270 
0C230 
0C290 
00300 
0C31C 
00320 
00330 
003^0 
CC350 
00360 
C0370 
0038C 
0C390 
OO^OC 
00410 
00420 
00430 
00*i40 
00450 
00460 
00470 
00480 
0049C 
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(2nd  card) 

TA7 
.U6 

no.  40 

690S1 

Bleick 

Orbital  transfer  with 
minimum  fuel . 

IS 

DEC  8f          3  05  7  JJ 

TA7  69081 

.U6    Bleick 

no. 40      Orbital  transfer  with 
minimum  fuel. 


genTA      7  U6  no.  40 

Orbital  transfer  with  minimum  fuel  / 


3  2768  001  61500  8 
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